Elliptic Rydberg States as Direction Indicators 
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The orientation in space of a Cartesian coordinate system can be indicated by the two vectorial 
constants of motion of a classical Keplerian orbit: the angular momentum and the Laplace-Runge- 
Lenz vector. In quantum mechanics, the states of a hydrogen atom that mimic classical elliptic 
orbits are the coherent states of the SO(4) rotation group. It is known how to produce these states 
experimentally. They have minimal dispersions of the two conserved vectors and can be used as 
direction indicators. We compare the fidelity of this transmission method with that of the idealized 
optimal method. 
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I. UNSPEAKABLE QUANTUM INFORMATION 

Information theory usually deals with the transmis- 
sion of a sequence of discrete symbols, such as and 1. 
Even if the information to be transmitted is of continu- 
ous nature, such as the position of a particle, it can be 
represented with arbitrary accuracy by a string of bits. 
However, there are situations where information cannot 
be encoded in such a way. For example, the emitter (con- 
ventionally called Alice) wants to indicate to the receiver 
(Bob) a direction in space. If they have a common co- 
ordinate system to which they can refer, or if they can 
create one by observing distant fixed stars, Alice simply 
communicates to Bob the components of a unit vector n 
along that direction, or its spherical coordinates and 
4>. But if no common coordinate system has been estab- 
lished, all she can do is to send a real physical object, 
such as a gyroscope, whose orientation is deemed stable. 

In the quantum world, the role of the gyroscope is 
played by a system with large spin. For example, Al- 
ice can send angular momentum eigenstates satisfying 
n • J|-0) = j\ip). This is essentially the solution proposed 
by Massar and Popescu [l| who took N parallel spins, po- 
larized along n. This, however, is not the most efficient 
procedure: for two spins, a higher accuracy is achieved 
by preparing them with opposite polarizations 0j. For 
more than two spins, optimal results are obtained with 
entangled states y, Q • 

The above discussion can be generalized to the trans- 
mission of a Cartesian frame. If N spins are available, 
one can encode a Cartesian frame in an entangled state 
of these spins, as in However, a more accurate trans- 
mission is then obtained if Alice uses half of the spins 
to indicate the s-axis, and the other half for her y-axis 
0. In this case the x and y directions found by Bob 
may not be exactly perpendicular and some adjustment 
will be needed to obtain Bob's best estimate of the x 
and y axes. Finally, the z-axis can be inferred from the 
estimates of the x and y axes. 

However, it is not possible to proceed in this way if 
a single quantum messenger is available. The optimal 



transmission of a Cartesian frame by a hydrogen atom 
(formally, a spinless particle in a Coloumb potential) was 
derived by Peres and Scudo • The results of [j| can also 
be used for a hydrogen atom if one considers the angular 
momentum eigenstates to be those of the atom, rather 
than those of N spins. 

In this paper we show how to transmit a Cartesian 
frame by using the elliptic Rydberg states of a hydrogen 
atom. These are the quantum mechanical analog of a 
classical Keplerian orbit, and it is known how to produce 
these states experimentally. Elliptic Rydberg states, just 
as their classical counterparts, define three orthogonal 
directions in space, and thus are natural candidates for 
encoding a Cartesian frame. 

In the following section we discuss the properties of 
quantum elliptic states. Section III deals with the trans- 
mission of one direction by means of them, and in Sec. IV 
we use them to transmit two orthogonal axes (and thus 
a Cartesian frame). In Sec. Ill and IV the detection pro- 
cedure is based on S0(3) coherent states as in 7]. S0(4) 
coherent states are employed in Sec. V to produce a posi- 
tive operator valued measure (POVM) which enables the 
use of elliptic states for the transmission of two direc- 
tions that are not orthogonal. Even when these states 
are used by Alice to transmit two orthogonal axes, the 
two directions found by Bob are not necessarily perpen- 
dicular and further adjustment is needed, as explained 
above. As shown in the Appendix, these adjustments in- 
crease the fidelity of transmission of two orthogonal axes. 
However, a higher fidelity is achieved with a POVM based 
on the S0(3) rotation group, especially when the energy 
quantum number n is large. 



II. CONSTRUCTION OF AN ELLIPTIC STATE 

A classical bounded Keplerian orbit in a potential 
— k/r can be defined by its constants of motion: the en- 
ergy E < 0, the angular momentum L which is an axial 
vector perpendicular to the plane of the orbit, and the 
Laplace-Runge-Lenz (LRL) vector 



K = {-2H)- 1/2 (p x L - fikr/r) , 



(1) 
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where fx is the particle's reduced mass, and we introduced 
a prefactor (— 2H)~ 1 / 2 for later convenience. This pref- 
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actor, which is a constant of motion, does not appear in 
the usual definition of the LRL vector The Hamil- 
tonian is H = p 2 /2m — k/r. We consider only bounded 
motion for which the energy E, which is the numerical 
value of H, is negative. The classical orbit is then an el- 
lipse, and the LRL vector is a polar vector directed along 
its major axis. It satisfies: 

L • K = 0, (2) 

and 

K 2 + ^L 2 = -i?k 2 l2E. (3) 

Because of these relations, only five out of the seven con- 
stants are independent. They uniquely determine the 
shape and orientation of the ellipse. Its eccentricity is Q 



e = \K\^/^2E/nk. (4) 

We now turn to the quantum version. We use natural 
units, fj, = k = h = 1, so that the energy levels for 
bound states are E = — l/2n 2 , where n is an integer. 
The operator K is defined by 

K=(-2F)- 1 / 2 (i(pxL-Lxp)-r/r). (5) 

Note that H commutes with (i(px L — L x p) — r/r). 
The commutation relations for the operators L and K 
are 

[Li,Kj] = ie ljk K k , (6) 

[Ki,Kj] = ie ljk L k . (7) 

Together with [L^, Lj] — ieij k L k , these are the commu- 
tation rules of infinitesimal rotations in four dimensional 
Euclidean space, which leave the n-th energy level sub- 
space invariant [jjj]. The coherent states of SO (4), i.e. 
the states for which the dispersion of L 2 + K 2 is mini- 
mal, can be built from the coherent states of SO (3), since 
SO(4) =SO(3)xSO(3). 

Define 

Ji = i(L-K) and J 2 = ±(L + K). (8) 

These two operators have the commutation relations of 
two independent three-dimensional angular momenta: 

[Jii, Jxj] = ieijk Jik, (9) 



[Jaii J"2j] = itijkJsk, (10) 



Instead of the classical equations (0) and Q we now have 

m 

L K = K L = 0, (12) 

and 

L 2 + K 2 = -1 - 1/2H = n 2 - 1, (13) 

where the last form of the equality holds for energy eigen- 
states. In the classical limit, n ^> 1, Eq. Ijl3|l reduces to 
the classical one ©. 

In the rest of this article we consider only energy eigen- 
states. Owing to Eq. (jHJ, we have 

3i(h+l)=h(h + l) = l(n 2 -l). (14) 

where j\ and j% are the quantum numbers referring to 
the operators J 2 and J 2 respectively. It follows that j\ 
and j2 are equal: j\ = ji = j, and j(j + 1) = j(n 2 — 1), 
so that 

j = I(n-l). (15) 

The coherent states of a three-dimensional angular mo- 
mentum will be denoted by | J, u). They obey u- J| J, u) = 
j | J, u) for an arbitrary classical unit vector u. For the 
coherent states the dispersion AJ = ((J 2 ) — (J) 2 ) 1 / 2 
is minimal: (AJ) 2 = j. In particular, AJ U = 0, and 
AJi = \J j /2, where J u = J • u and Ji = J • v with 
v _L u. The coherent states of SO (3) are obtained by a 
rotation of a fiducial coherent state | J, z) , 

\J,u e4> ) =e-* L ^e- iL y e \J,z). (16) 

The coherent states of SO (4) are now obtained as di- 
rect products of coherent states for each of the SO (3) 
subgroups, 

1 71, UiU 2 > = |Jl,Ui) <g) |J 2 ,u 2 ). (17) 

where the unit vectors Ui and u 2 are again classical. The 
coherent state |Ji,Ui) obeys 

ui • Ji|Ji,Ui) = j|Ji: u i> = \{ n - l)l<A> u i), (18) 

and likewise we have 

u 2 • J 2 |J 2 ,u 2 ) =j\J 2 ,u 2 ) = i(n-l)|J 2 ,u 2 ). (19) 

As from now, we shall omit the symbols n, J\, and J 2 in 
state vectors, since the quantum numbers 7£,ji,j 2 have 
fixed values, related by Eq. 114|) . For example the state 
\nlm) which obeys H\nlm) = n\nlm), ~L 2 \nlm) = 1(1 + 
l)\nlm), and L z \nlm) = m\nlm) will be written simply 
as \lm), \ J, u) becomes |u), and |J, z) becomes \jj), etc. 
The symbol j will always denote the fixed value j = 
i(n-l). 

Owing to Eq. JSJl, the dispersion of L 2 + K 2 is minimal 
for coherent states 1 1 11 : 



[Ju, J 2 j] — 0. 



(11) 



(AL) 2 + (AK) 2 - 2[(AJx) 2 + (AJ 2 ) 2 ] = 2(n - 1). (20) 
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To obtain the expansion of the coherent state |n, uiu 2 ) We then use the angular momentum addition formula, 
in the familiar nlm basis, we first expand each of the 
in Eq. flJJ: 



|ui) = £ D L(9i4>i)\ 



jm) 



1 = 1,2, (21) 



m=-j 



where the -0^(6^^) are related to the usual rotation 
matrices LL2I: 



Dl 



v^(cf>eo) mj , 



(22) 



27 



1/2 



+ rn 

Substitution into Eq. I|17|) gives 



(cosf) J+m (sin|y- m e- mi0 . (23) 



2j t 



\jrrn) ® |jma) = E E J Z ™>> ( 25 ) 



i=0 ra=-/ 



|uiu 2 ) = E E D mi(^^) D m 2 {824>2)\jm 1 )<g)\jm2). where C^f 2OT is the Clebsch-Gordan coefficient [3 

which vanishes for m 7^ mi + TO2. Combining Eqs. I|24|) 
(24) and JSKJ, we finally get 



mi=-j m 2 = -J 



2j I _ 

l U l U 2)=EE I E E D LS e ^) D L{ e ^) CJ m,rn^> I I'"')- 
Z=0 m— — l \mi=—j rri2=—j 



r 



(26) 



The classical orbit that corresponds to the coherent state 
|uiu 2 ), in the limit of large n, can be obtained as follows. 
From (JSJ, we have 

L=(Ji+J 2 ) and K = (J 2 -Ji). (27) 

Let C be half the angle between ui and u 2 , i.e. sin£ = 
I Hi x u 2 |/|ui + U2I, and define three orthogonal classical 
unit vectors 



1 = 



Ui + u 2 
|ui + u 2 | 



U2 Ul 

|u 2 - Ui| 



(28) 



and w = i x k. Denoting by an arbitrary vector 
orthogonal to Ui, we have 

Ji = (Ji •Ui)ui + (Ji ■ui_ L )uix, (29) 



we get 



(ui|(u 2 ■ Ji)|ui 



(u 2 |(ui-J 2 )|u 2 ), (33) 
j ui • u 2 = j cos 2(. (34) 



Noting that 



and 



|m + u 2 | = 2| cos C|, 



|ui - u 2 | = 2|sinC|, 



(35) 



(36) 



we obtain from (|34|l the expectation values of the com- 
ponents of K and L along the directions of k, £ and w, 
for the coherent state |uiu 2 ) 



(K k ) = (n-l)sinC, 



(37) 



Jl ■ U 2 = (Ji ■ Ui)(Ui • U 2 ) + (Ji • Uix)(ui_L • u 2 ). (30) 

Then from 

(ui|(ui- Ji)|ui) =j, (31) 

and 

(ui|uij.-Ji|ui)=0, (32) 



{Li) = (n - l)cosC, 



(38) 



where Kk = k • K, etc. In the perpendicular directions 
the expectation values vanish: 



(Kt) = (K w ) = (L k ) = (L w ) = 0. 



(39) 



From Eqs. (|37|I - H39(I we see that in the limit of large 
n, the coherent state |uiu 2 ) corresponds to a classical 
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elliptic trajectory in the kw plane with the LRL vector in 
the k direction, the angular momentum in the £ direction, 
and eccentricity e = (Kk)/(n — 1) = sin£, which is the 
quantum mechanical analog to Eq. Q. The unit vector 
w is parallel to the minor axis of the ellipse. 



for each I. Since \B) is a direct sum of vectors, one for 
each value of I, then likewise U (ipOfi) is a direct sum with 
one term for each irreducible representation: 



(47) 



III. TRANSMISSION OF ONE DIRECTION 

We now turn to the use of elliptic wave functions as 
direction indicators. Consider two observers (Alice and 
Bob) who do not have a common reference frame. Alice 
wants to indicate to Bob her z-axis by using a hydrogen 
atom in a Rydberg state. In the next section, we shall 
likewise discuss the transmission of two orthogonal axes 
by a single hydrogen atom. We use as much as possible 
the same notations in both sections. Alice's signal is 



2j / 



1=0 m=-l 



} J ai m \lm), 



where 



2j 



E E 



O/r, 



(40) 



(41) 



i=o 



Bob's detectors have labels ip9cj) which indicate the un- 
known Euler angles relating his Cartesian axes to those 
used by Alice. The mathematical representation of his 
apparatus is a POVM 



dE(ip0<j>) = 1, 



(42) 



where 



dE(ipecf>) = d^u{^ecj>)\B){B\u\^e4>), (43) 



and d^otp = sin 9d^d9dcj}/8Tr is the SO(3) Haar measure 
for Euler angles [l^- As usual, U(ip9cf>) is the unitary 
operator for a rotation by Euler angles ip9<fi, and \B) is 
Bob's fiducial vector defined as in 



2/ 



\B) =]T\/27TT J2 b im\lr 



1=0 



where for each I, 



E fa 



(44) 



(45) 



Note that Eq. (|40|l was written with Alice's notation, 
while Eq. (|44|l is in Bob's notation (recall that they use 
different coordinate systems). 

Optimizing the transmission fidelity, defined by 
Eq. gSJl below, leads to 



hm = aim ^ E l a '"| 2 ^ 



-1/2 



(46) 



where the V ^ l '{iji6^) are the usual irreducible unitary 
rotation matrices |12(. A generalization of Schur's lemma 
14] confirms that Eq. (|43|l is indeed satisfied, owing to 
the coefficients \J21 + 1 in Eq. 

The fidelity of the transmission of a single direction is 
defined as usual: 



F = (cos 2 (w/2)) = |(1 + (cosw)), 



(48) 



where u> is the angle between the direction indicated by 
Alice and the one that is estimated by Bob. If Alice 
indicates her z-axis, we thus want to maximize (cosa^). 
Following the method of Peres and Scudo we define 
Euler angles a[3j whose effect is rotating Bob's Cartesian 
frame into his estimate of Alice's frame, and then rotating 
back the result by the true angles from Alice's to Bob's 
frame. The angles a(3j thus indicate Bob's measurement 
error. Since in this case Bob's estimate refers to Alice's 
z-axis only, the angle to z is identical to the second Euler 
angle (3: 



(costo z )= J d a01 \(A\U(af3~f)\B)\ 2 cos 13. (49) 

Let us examine two extreme cases. First, we take Al- 
ice's vector to be a circular state, with null eccentricity 
(sin C = 0), i.e., 



with 



\A) = \U), 



I = 2j = n - 1. 



(50) 



(51) 



Bob's vector is obtained from (|46|l . which in this case 
gives 



\B) = V2n-1 \ll). 



(52) 



We then have 1# 



(A\U(apy)\B) = V2^~T e^ n -^ a+ ^ cos 2(n ~ 1] ^ (53) 
Inserting the last equation into Eq. I|49(l gives 

(coBU z ) = (n-h) ( sin/3d/3cos 4( "- 1) (/3/2)cos/3, (54) 



(55) 



= (n- I) /n. 



The "infidelity" 1 — F, whose typical meaning is Bob's 
mean square error Q , is 



i(l - (coslu z )) = l/2n. 



(56) 
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This result is identical to the one obtained by Massar 
and Popescu Q , when the number of their parallel spins 
is taken to be N = 21 = 2n — 2. 

The other extreme case corresponds to a classical el- 
liptical orbit with unit eccentricity, so that L = 0. Let 
K lie in the z direction. The corresponding quantum 
state, denoted by \K,z), is an extreme Stark state with 
(L z ) = and maximal (K z ): 



\K,z) = |-z) 



|z) = \j, ~j) ® \jj). 



(57) 



This state satisfies L z \K,z) = and it is also an eigen- 
state of K z : 

K z \K,z) = (J 2z -Ju)\j,-j)®\jj), (58) 
= (n-l)\K,z), (59) 

owing to 

( - Ju)\j, -j) ® \jj) = -Ju \j, ~j) ® hz \jj) ■ (60) 
In the nlm basis we have 



\K,z) 



2j I 

EE 

;=o m=-l 



2j 

CU! m2 Jlrn)=J2 C -%\ 10 )' ( 61 ) 

1=0 



since mi = — j and m 2 — j. The fidelity of transmission 
by this state will be evaluated at the end of this Section. 

Both the circular state and the extreme Stark state 
are coherent states of SO(4), but only the circular state 
is also an angular momentum coherent state. Moreover, 
the circular state is symmetric, (ll\r\ll) = 0, while the 
extreme Stark state state is not. This can be seen from 



(nlm\r\nlm) — ^(nlm\K\nlm) . (62) 



Let us examine which one of these states gives better 
results when used by Alice to transmit the directions of 
her z-axis. The overlap between two angular momentum 
coherent states is LL5| 



|(u 1 |u 2 )|^=cos^( X /2), 



(63) 



where x is the angle between the directions of Ui and u 2 . 
It is noteworthy that the overlap between two extreme 
Stark states is the same, as we will see shortly. First, a 
rotation of the \K,z) state by angles (9<fi) gives 



\K,z), 



(64) 



where again the operator e %L ^e %L y performs an ac- 
tive rotation of the vector \K, z). Using 127|) we have 



\K, ues) = e 



-iJ2z4>p — iJ2y8 



|z), (65) 
*>, (66) 



owing to Eqs. 10 and Q). Thus the rotated extreme 
Stark state is just 



\K, u fl 



(67) 



where the SO(3) coherent states \ug^) are defined as in 
Eq. H21fl . This Stark state is an eigenstate of u • K with 
the maximal eigenvalue n — 1, and it satisfies k = u, as 
can be seen from Eq. (|28|1 . The overlap between two such 
states, K^u'l^u")! 2 , is 



K- 



(68) 



which by using (|63|l is just: 

cos 4jl (x/2) cos 4j2 (x/2) = cos 4( "- 1) (x/2), (69) 

where x 1S the angle between the vectors u' and u". 

Such a simple expression cannot hold for the overlap 
of two generic elliptic states whose eccentricities are not 
or 1. Let a generic elliptic state 



UiU 2 ) 



u 2 , 



(70) 



be an elliptic state with eccentricity < e < 1. Unlike 
the e = 1 and e = cases, this state does not define one 
direction, but two independent ones Ui and u 2 . If it is 
rotated by Euler angles a[3j the result is 



-iL z a —iLy(3 —iLrf 



e-^lm u 2 ) = Ux\xli) 8 U 2 \u 2 ). (71) 



where 



_ p — iJizOi p -iJlyl3 -iJizl 



(72) 



and likewise for C/ 2 . To obtain this result we have 
used l|27|) and the commutation relations © and Q) . The 
rotation e~ ji * a e~ lLy ' 3 e -lLz7 opens an angle %i between 
the classical vectors Ui and i?(a/?7)ui, and an angle xi 
(which is generally different from xi) between u 2 and 
i?(a/37)u9. Here R(a/3j) denotes the classical rotation 
matrix 81. It follows that 



\ ( Xi X2 
UiU 2 ) = {cos — cos — 



2(n-l) 



(73) 



Generally, both xi an d X2 are different from the angle 
between the directions k and k' = i?(a/37)k, or between 
the directions i and t! = R(a/3j)£. 

We now calculate the transmission fidelity for the case 
where Alice sends an extreme Stark state \K, z). Since 
| A) contains only m = terms, so does Bob's fiducial 
vector 



Wc thus have 



bim — a;o(|a;o| 2 ) 1//2 <5mO- 



blm — 8m0 ( a io/| a Zo|), 



\B) = Y,V2l + l (o w /|oio|)|/0). 



(74) 
(75) 

(76) 
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In order to determine (cosu; z ) in Eq. I|49[) . we note that 

{A\U(a0i)\B) = J2 V^l+Ta^btoilOlV^ia^lO), 
1=0 

n-l 

= ^^/2^Tl|a ^ o|d ( ^ ) ) ( ^ 9). (77) 



2=0 



We insert this expression into i|49|). The result, obtained 
by using Eqs. (19)-(21) of ref. U is 



(C0SW Z ) = |ajO Ctjfeol, 

kl 



(78) 



where is a real symmetric matrix whose non- 

vanishing elements are 



(79) 



and a;o = ^-jjo- The results are summarized in Fig. 1, 
in which the mean square error is plotted versus n. The 
| if, z) state gives better fidelity than the circular state 
but for n > 3 its fidelity is substantially less than 
optimal 0)0 and goes asymptotically to l/(4n— 2). This 
raises the question whether it is possible to build a "nat- 
ural" POVM by setting Bob's vector to \B) = VN\K, z), 
so that POVM elements are 



N\K, u e0 )(if, u e 



(80) 



where |if, ug^) was defined in Eq. 1)66(1 and N is a nor- 
malization factor. Unfortunately, |if, z) contains a su- 
perposition of all values of I, as can be seen from 1)6 
Thus | if, z) does not belong to one irreducible subspace 
of the representation of the SO(3) rotation group. As a 
result, the operator 



B= I d e 0|if,U00)(if,ug 



(81) 



is not proportional to the identity, but is a block-diagonal 
matrix with different blocks for each irreducible repre- 
sentation of the rotation group. Moreover, the resulting 
POVM includes an element which corresponds to the ab- 
sence of any answer, thus reducing fidelity. A natural 
POVM which uses the SO(4) group will be discussed in 
Sec. IV. 

The direction of the minor axis of a classical nondegen- 
erate ellipse is that of L x K. A quantum ellipse also has 
this property. Taking Alice's state as a quantum ellipse 
with eccentricity < e < 1, with both k and £ lying in 
the xy plane so that w = z, the resulting fidelity can be 
compared with the cases where k or t points along the 
z-axis and the eccentricity of the ellipse is or 1, respec- 
tively. The fidelity for transmission using the semi-minor 
axis reaches a maximum at eccentricity of about e = 0.7 
(a different eccentricity for each value of n) . A compar- 
ison of the mean square error for using the three options 
is given in Table 1. 



A. Comparison between elliptic wave functions and 
optimal wave functions 

We shall now compare the extreme Stark state with 
Alice's optimal vector for the transmission of one axis as 
calculated in 0. They are both eigenstates of L z with 
m = 0, and since Eq. (|78|l holds, we will present them in 
the notation: (|a 00 |, Ko|, |a 2 o|, kn-i,o|), where a ;o = 
C-jjO as before. For n = 3 we have 



while Alice's optimal state is 



(82) 



(83) 



Thus for n — 3 the overlap between the extreme Stark 
state and the optimal state is 



\(K,z\A 



opt) 



0.993491. 



(84) 



Both states give almost the same fidelity for transmission 
of one axis. For higher values of n, they become more 
and more different. For n = 10 the overlap is 



\(K,z\A 



opt/ 



0.76406. 



(85) 



The various components are given in Table 2. We see 
that the extreme Stark state has coefficients peaked at 
lower values of I than the optimal state. 



IV. 



TRANSMISSION OF TWO AXES 



Alice now wants to transmit a Cartesian frame by in- 
dicating the directions of two axes, the third one being 
inferred from them. Which elliptic state is optimal? Ob- 
viously, states with e = and e = 1 will not do in this 
case, since they define only one direction. We have to 
find the optimal eccentricity. Let 



Aif ± = V(Aif^) 2 + (AK w f 



AL ± = ^{AL k y + {AL V 



(86) 



(87) 



where 



AK f 



owing to Eq. (|39|l . We define similar expressions for the 
other components. When we want to transmit k, namely 
the direction of the classical LRL vector, then a smaller 
Aifj_/(iffc) improves the fidelity. A similar argument 
holds for the transmission of t. Thus when transmitting 
two axes, a heuristic guideline is to look for states that 
satisfy 



Aif I 



ALj_ 

WW' 



(89) 
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A straightforward calculation gives 

AK W = AL W = ^/i(n-l), (90) 

A^ = yi(^l)sinC, (91) 

AL fc -^/i(n-l)cosC (92) 

Together with Eqs. i|37[l and l|38|) . this gives an equation 
for the eccentricity, 

Vl + sin 2 C ^ y/1 + cos 2 C (Q3) 
^/2(n- IJsinC ~ \/2(n - IJcosC' 

Therefore we expect that the optimal eccentricity is ap- 
proximately 

e = sin C = cos C = 1 / V2. (94) 

More accurate numerical results are given below. 

We now evaluate the fidelity for the transmission of two 
axes. Alice uses an elliptic state with k = x and £ = y 
(the unit vectors in the x and y directions respectively). 
The eccentricity e = sin £ has to be optimized. Recall 
the £ is defined to be half the angle between Ui and 
u 2 . The definitions of k and £ are given in Eq. Q28JI. 
Thus in order to meet the above requirements we set in 
| A) = \ug 1< p 1 ,ug 2 ^ 2 ) the parameters 61 = 62 = tt/2, and 

01 = |tt - C, 02 = |tt-C- (95) 

Fidelities now must be defined for each one of the axes. 
Note that coswfc (for the kth axis) is given by the corre- 
sponding diagonal element of the orthogonal (classical) 
rotation matrix. For the transmission of the x and y 
axes, we thus need @ 

(cos uj x + cos ujy) = ((1 + cos /?) (cos(a + 7)). (96) 

We expand |A) and \B) as in (gUJ) and gU- Bob's 
optimal fiducial vector is still given by (|4^|) , and (cos u>x + 
cos ujy) is calculated using equations (23)-(26) of ref. Q. 
The mean square error per axis is plotted in Fig. 3 for 
n = 5, 10, and 20. The error is minimal at e s» 0.708 for 
n = 5, at e w 0.704 for n = 10, and e « 0.674 for n = 20. 

Note that the shape of the curve flattens with increas- 
ing n, so that the minimum is hard to find numerically. 
The intuitive explanation is that in the limit of large n, 
as if Alice were to transmit a "classical atom," i.e., a 
classical two body Kepler system, then the direction of 
the classical angular momentum and LRL vectors could 
be found irrespective of the eccentricity. Therefore, the 
transmission accuracy would be the same for any eccen- 
tricity that is not close to zero or one. 

The deviation of the optimum from e = l/v2 was 
expected, since transmission of the k direction (e = 1) 



achieved higher fidelity than the transmission of the £ 
direction (e = 0). Thus the ellipse with optimal ec- 
centricity for transmission of two axes is biased to give 
AL±/L < AKx/K in order to compensate the differ- 
ence and make the contribution to the error from the k 
direction about equal to that from the £ direction. 

Elliptic states give results very close to the optimal 
ones. The mean square error for transmission of two axes 
by elliptic states with optimal eccentricity is compared to 
the optimal results in Table III. 

V. POVM FOR SO(4) 

We now construct a POVM based on the SO (4) group 
and use it in order to transmit two axes. This POVM is 
naturally built with the SO (4) coherent states which are, 
as we have seen, direct products of two SO (3) coherent 
states. We shall use for each one of the SO(3) subspaces 
the notation 

|#0) u = V27TT f/(W)|u), (97) 

and 

dE(ip6ct>) = d W |W) u (W|u, (98) 

where u labels the direction to be transmitted, and 
d^e<t> = s'mdddd(j)dip/8ir 2 as in Eq. (|43|l . By applying 
Schur's lemma to each of the SO (3) subspaces we have 

J I dEx^e^t) ® dE 2 (ij 2 6 2 (j )2 ) = li ® 1 2 = 1. (99) 

We are now ready to discuss the transmission of Alice's 
x and y axes by means of an elliptic state. We take 

I A) = |xy) = |x) ® |y). (100) 

This equation was written in Alice's notation. We also 
define a fiducial vector for Bob 

|S) = (2j + l)|x)®|y), (101) 

written in Bob's notations. Thus the POVM element is 
constructed from the vector 

I^kMx ® IV'2^20 2 ) y = (2j + l)(t/i ® U 2 )\B). (102) 

The result of Bob's measurement consists of two sets 
of Euler angles, V'i^i'^i an( i 4 ! 26 2 4> 2 . The first one gives 
Bob's estimate of the active rotation needed to bring his 
x-axis to Alice's x-axis. Likewise, the second set gives 
Bob's estimate of the active rotation needed to bring 
his y-axis to Alice's y-axis. The detection probability 
of these sets of angles is 

dP(V>i...«fe) = dfro^dfrihto I (8 U 2 \B)\ 2 . (103) 

Recall that Eq. (|100|) was written in Alice's notations, 
while (|101fl is in Bob's notations. To compute the result 
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explicitly, we need a uniform system of notations. For 
this we introduce, as in 0, the Euler angles £r]( that ro- 
tate Bob's xyz axes into Alice's axes. (The Euler angle 
£ should not be confused with the eccentricity parame- 
ter introduced in Sec. II.) The unitary operator U (£?K) 
represents an active transformation of Bob's state vec- 
tors to the corresponding state vectors of Alice's system. 
Therefore, U(£r)£) is also the passive transformation from 
Alice's notations to Bob's notations. Written in Bob's 
notations, Alice's vector \A) becomes U(£r)£)\A) so that 
in Eq. (A\ becomes (A\U(£r)C)i. Owing to the 

commutation relations © and Q 

UfcvQ = e- iL ^e~ iL ^e' iL ^, (104) 

= tfi(£0®£M£0, (105) 

where again U\ and U 2 are defined as in Eq. J75J. Wc 
thus have 

t/(^C)|u lU2 ) = £M£OK) ® C/ 2 (^C)|u 2 ). (106) 
Let us therefore define 

Vx{a x frti) = uUfriQ E/i(V>i (107) 

and 

U 2 (a 2 (3 2l2 ) = uUtvQ U 2 (yj 2 6 2 fo). (108) 

We shall henceforth use the left hand sides of Eqs. 1)107(1 
and l|108|) as the new definitions of the symbols U\ and 
U 2 . As before, the Euler angles ct\fi\^\ have the effect of 
rotating Bob's ir-axis into his estimate of Alice's x-axis 
and then rotating back the result by the true rotation 
from Alice's to Bob's frame. The action of the Euler 
angles a 2 (3 2 j 2 is similar for the y-axis. Thus the Euler 
angles otifii^i indicate Bob's measurement error, and the 
probability of that error is 

dP(ai...j 2 ) = d ai ^ 1 d a2 ^ 2 \(A\Ui^U 2 \B)\ 2 . (109) 

Note the similarity with Eq. 1(103(1 . The difference is that 
((103(1 referred to the probability of detection of a partic- 
ular set of Euler angles, while l|109|) gives the probability 
of error in that detection. 

The transmission mean square error per axis is, as in 
Eq. 



where we have used (|109ll and Schur's lemma for the 
second set of angles, namely 

(2j + 1) J d a2P2l2 U 2 \y)(y\Ul = 1 2 . (112) 

The evaluation of Eq. 1(111(1 is identical to the one per- 
formed in Eq. I(54|) . with n replaced by h{n + 1) every- 
where, and we get 

(coslu x ) = (n-l)/(n + l). (113) 

Likewise 

(cosujy) = (n- l)/(n + 1). (114) 

Thus the infidelity (mean square error) per axis is 

|(1 - (cos^)) + i(l - (cosw v )) = l/(n + 1) (115) 

In ref. it was found that the optimal POVM (not re- 
stricted to elliptic states) for transmission of two axes 
using a hydrogen atom is of the form given by Eq. 1(43(1 . 
It was shown that using this POVM, the infidelity per 
axis for Alice's optimal signal approaches l/(3n) asymp- 
totically. Using SO(4) instead of SO(3) as in Q, we ob- 
tain an infidelity per axis that is exactly l/(n + 1) for all 
values of n. As shown in the Appendix, an adjustment 
procedure to obtain to orthogonal axes will further de- 
crease the mean square error by a factor which for large 
values of n tends to 3/4. 

The SO(4) POVM also enables the transmission of two 
directions which are not orthogonal, by means of a sin- 
gle hydrogen atom in an elliptic state. To transmit the 
directions of two general unit vectors vi and v 2 , Alice's 
prepares the elliptic state 

|A) = |viv 2 > = |vi)®|v 2 >, (116) 

(in her notations) while Bob's vector is (in his notations) 

|fl) = (2i + l)|vi)®|v 2 >. (117) 

As before, the infidelity for each direction is l/(n+ 1). It 
should be noted that transmission of two non-orthogonal 
directions with one hydrogen atom is not possible with 
the SO(3) POVM. 



R = j(l - COSWa) + |(1 - COS£J„), (110) 

where u x and uj y are the angles between the true and es- 
timated directions of the s-axis and y-axis, respectively. 
Since Bob infers the direction of the ir-axis from the an- 
gles ipi6i4>i, the value of cosa^ depends only on a 1 /3 1 7 1 . 
Likewise, the value of cosw y depends only on the angles 
a 2 /3 2 j 2 . We have 

(cosuj x ) = / d Ql/3l7l |(x|J7i|x)| 2 cos^, (111) 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We have shown how elliptic Rydberg states can trans- 
fer information on the orientation of one direction, or 
more generally that of a Cartesian frame. For increasing 
values of n, the fidelity obtained for a single direction falls 
rapidly below the optimal ones. However, for a Carte- 
sian frame the results are very close to the optimal ones. 
Furthermore, elliptic states have the advantage of being 
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experimentally accessible, while preparation of the opti- 
mal states seems much more difficult. Note that we have 
assumed Alice and Bob have the same chirality. If their 
chiralities are opposite, then when angular momenta are 
used for the transmission, the direction inferred by Bob 
should be reversed (because directions are polar vectors 
while angular momentum is an axial vector). However, 
the LRL vector is also a polar vector, thus even if Bob 
and Alice have opposite chiralities, the direction inferred 
by Bob is correct. We have also shown how elliptic Ry- 
dberg states can be prepared to encode two arbitrary 
directions, when the measurement is based on the SO (4) 
rotation group. 

Work by AP was supported by the Gerard Swope Fund 
and the Fund for Promotion of Research. Work by NHL 
was supported by a grant from the Tcchnion Graduate 
School. 



t>2 = 02 - 9 71 "- 



and 



cosw y = sin 02 cos 02 ~ 1 — \Q\ — \<j>2- 



(A7) 



The deviation angles are given by cosui x = f x ■ x and 
cos a; y = f y ■ y, namely 



cosw x — sin 6>i cos 0i « 1 — — \(j>i, (A8) 



(A9) 



Let g denote the infidelity per axis before the adjust- 
ment. The infidelities for both axes are equal, thus 



g = |(1 - (coswz)), 



(A10) 



\ I ' {Ol + c^dpi = \{el + ^), (All) 



APPENDIX A: REDUCTION OF ERRORS BY 
ORTHOGONALIZATION 

As we have seen in Sec. V, Bob's estimates of Alice's 
x and y axes may not be exactly orthogonal. The prob- 
ability for the estimate of the x axis to have an angular 
error u x , as can be seen from Eq. i|t)3|) . is 



p{uj x ) oc cos^™ \uj x /2) 



(Al) 



and likewise for the y axis. Thus for large values of n, the 
error probability distribution will be highly peaked. We 
now calculate the gain in fidelity achieved if Bob performs 
a simple orthogonalization of his two estimates f x and f y , 
by rotating the two vectors in their plane by the same 
angle, so that they become orthogonal. 

Let us define two pairs of spherical angles that give the 
position of the estimated directions with respect to the 
(unknown) true axes. These positions are given by 



r x = (sin Q\ cos 0i , sin 9\ sin 0i , cos Q\ ) , ( A2) 



and 



where 



fullfills 



dp l = p(9i, 4>i) sin 9i dOi dfa 



dpi = 1. 



(A12) 



(A13) 



Equivalently, we can write the infidelity in terms of 62 
and 4>2 as 



dp % = 



4\ u 2 



(A14) 



In first order we have, by combining (|A2|I and (|A3() with 
the definitions ((A"6|l and (|A7|) . 



(1,^1, 



(-02,1,-02), (A15) 



and the angle Q between them is given by 



cos ft = f x ■ r y 



»2- 



(A16) 



(sin 62 cos 02 , sin 62 sin 02 , cos 62 ) ■ 



(A3) 



The probability distributions will be denoted Pi(0i,<fii). 
In the limit of large n, the deviation angles uj x and uj v 
are small. Hence the distribution are centered as 



Px 



= p{9 X - ±7T, 0l), 



Py = /?(#2 - ^71", 02 - ±71"), 



(A4) 



(A5) 



where p(£,p) is peaked around (0,0). Here we used the 
fact that the SO(4) POVM gives probabilities of error for 
each axis which are identical and independent. Define 
new variables 



The bisector of f x and f y is given by the unit vector 
b = (fi +f 2 )/|?i + f-2 1 - Using l|A15|) and keeping only 
first order terms we have 



+ 2 ), 1 + i(0i + 2 ), -0i - h\ /V2, 

(A17) 



1- 2O1 ' ru,- 1 2 



where we used 

I? 1 + T2 



|fl+^|«^(l+|01-|0 2 ). 



(A18) 



We can also express the bisector b in terms of its spherical 
angles which we shall denote by (r, i^>) . Since the errors 
are small, we have ip « jTT, and it is convenient to define 



fa 



(A6) 



(A19) 
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Comparison of the two expressions for b gives 



i = 1^+^1^+92), £=i(^ + 4> 2 ). (A20) 

In first order, as Eq. i|A16(l shows, the orthogonaliza- 
tion consists in changing the angles 4>i irrespective of 9i, 
without changing the #j themselves. Hence, in first order, 
the procedure defines 



t>'i = <p- JTT, 



t>2 = V + \^-, 



i.e., 



b'l = 02 



= |(01 +02), 



(A21) 



(A22) 



where again <^> 2 = <j/ 2 — \k. The change in 6i is of higher 



2 

order, 6[ = 6i + 0(6 2 ,4> 2 ). The new infidelity per axis 



j"™ is 



new _ 1/^2 , 22\ _ 1(02 2 +6*1) 



3 ncw = |(0i +9i) 



(A23) 



Returning to (|A11|1 and (|A14|i . consider the integrals 
over <bj. Define 



P4,(4>) = / p{6,<j>)sm6d6 



(A24) 



Keeping in mind that the distributions for f>\ and </> 2 are 
identical, the 0-part of the infidelity per axis before the 
adjustment is 



W\) = -M) 



(A25) 



0?p <A (0i)p^(<^ 2 )d0i# 2 /4, (A26) 



,)d0i# 2 /4. (A27) 



The 0-parts of the infidelities after the adjustment, 
denoted by g^ ew , are 

9T W = Wi) = i(02 2 > = + 2 0i02 + 02>- (A28) 

The functions /5^(0i) and P0(0 2 ) are even, because the 
probability distribution p depends only on the angles u x 
or w y , which are independent of the sign of (f>i and 4>2 ■ 
Thus 

(0102) = / 0102 P (0l) P (0 2 ) #1 #2 = 0. (A29) 



With the help of Eq. (|A25|) we obtain 



9f w = Wi + &) = l M)- ( A3 °) 



Thus the (/>-parts of the infidelity are halved, 



new 1 / * ni\ 

g<f, = 2 50- (A31) 



As already stated, the angles 0i are unchanged in first 
order. Since the probability function p a (0i,^i) depends 
only on the angle lo x = r x - x, it is symmetric with respect 
to rotations around the x axis. A similar argument holds 
for the y-axis. Thus 



?1) = <0i> = {0l) = (02), 



and we have finally 



5 now = b- 



(A32) 



(A33) 
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n z = w z = I z = k 

5 e = 0.6963 e = e = 1 

?? = 0.193967 t? = 0.1 17 = 0.0573645 
10 e = 0.701261 e = e = 1 

= 0.0861934 i] = 0.05 r? = 0.0264067 

TABLE I: Eccentricities e and mean square errors jj for trans- 
mission of a single direction using z = w, z = £, or z = k, for 
n = 5 or 10. 



7i elliptic optimal 



5 


V 


= 0.14765 


'/ 


= 0.14465 


10 


V 


= 0.06822 


'/ 


= 0.06793 


20 


V 


= 0.03190 


'/ 


= 0.03088 



TABLE III: Mean square error 77 for transmission of two axes 
by an elliptic state with optimal eccentricity, and by the op- 
timal method Q for n = 5, 10, and 20. 



I 0123456789 

\K,z) 0.3162 0.4954 0.5222 0.4534 0.3365 0.2148 0.1167 0.0526 0.0186 0.0045 
Optimal 0.1825 0.3079 0.3767 0.4098 0.4130 0.3894 0.3422 0.2751 0.1923 0.0989 

TABLE II: Coefficients \aio\ for Alice's optimal state and for 
the extreme Stark state when n — 10. 
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Figure 1 

FIG. 1: Mean square error as a function of n for the transmis- 
sion of a single axis using the circular state (open circles), the 
extreme Stark state (squares), and the optimal state (closed 
circles). 




Figure 2 



FIG. 2: Mean square error (per axis) as a function of eccen- 
tricity for n=5, 10 and 20. 



